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ABSTRACT 

We investigate the dependence of cluster abundance n(> M,rci), i.e., the number 
density of clusters with mass larger than M within radius rd, on scale parameter 
rd- Using numerical simulations of clusters in the CDM cosmogonic theories, we 
notice that the abundance of rich clusters shows a simple scale invariance such that 
n[> {rci/rQ)"M,rci] = n(> M, tq), in which the scaling index a remains constant in a 
scale range where halo clustering is fully developed. The abundances of scale r^i clusters 
identified from IRAS are found basically to follow this scaling, and yield a ~ 0.5 in 
the range 1.5 < r^i < 4/i~^Mpc. The scaling gains further supports from independent 
measurements of the index a using samples of X-ray and gravitational lensing mass 
estimates. We find that all the results agree within error limit as: a ~ 0.5 — 0.7 in 
the range of 1.5 < rd < 4/i^^Mpc. These numbers are in good consistency with the 
predictions of OCDM {Qm = 0.3) and LCDM (Om + f^A = 1), while the standard 
CDM model has different behavior. The current result seems to favor models with a 
low mass density. 

Subject headings: cosmology: theory - galaxies: clusters: general - large-scale structure 
of universe 
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1. Introduction 

It is generally believed that the gravitational clustering of scale free initial perturbations 
is self-similar (e.g. Peebles 1980). The statistical descriptions of large scale structures, such as 
correlation functions and number densities, show somewhat of a scaling behavior. Among various 
observed objects, clusters probably arc the best for studying this gravitational scaling, because 
their distribution and clustering are dominated by gravitation and less "contaminated" by hydro 
processes. Therefore, observed clusters can directly be identified as massive halos of N-body 
simulation samples. A multi-scale identification of clusters will be useful to reveal the expected 
scaling. 

Traditionally, clusters arc identified from either observed or simulated distribution of galaxies 
or masses within a sphere of given radius ra, say, the frequently adopted Abell radius rd = 1.5h~^ 
Mpc where h = Hq/IOO km s~^ Mpc~^. However, a priori choice of rd is somewhat arbitrary, 
or dependent on observational limits. Clusters have already been identified and studied on very 
different radius r^i- For example, a wide range of radius covering from 0.01 - 4.3 Mpc has 
been used in the current determination of cluster masses (e.g. Wu & Fang 1996, 1997a and 
references therein). Therefore, a multi-scale identification of clusters is a necessary extension of 
the "standard" Abell radius rd = l.5h^^ Mpc. It will provide information, such as the scale 
invar iance of the mass functions of clusters, which can not be seen by one scale identification. This 
is the goal of this paper. 

As we know, cluster study has been significant for understanding the nature and evolution 
of cosmic structures on scales of ~ 1 - 10 Mpc. In particular, cosmological parameters can be 
well constrained by the observed properties of clusters, such as the abundances and correlations 
(Bahcall & Cen 1992 (BC); Jing et al. 1993; Jing & Fang 1994; Viana & Liddle 1996; Carlberg et al. 
1997a; Bahcall, Fan &; Cen 1997), the substructures (Jing et al. 1995), the luminosity-temperature 
relations (Oukbir Sz Blanchard 1997) and the gravitational lensing phenomena (Wu & Fang 
1996; Wu et al. 1997a). We believe that the multi-scale identification of clusters and the scaling 
invar iance would be able to add more constraints on relevant parameters. 

In §2, we summarize the theoretical background of the possible scaling of the multi-scale 
identified clusters. In §3, we present the result of identifying clusters on different scales in three 
popular structure formation models: 1) the standard cold dark matter model (SCDM), 2) low 
density, flat CDM model with a non-zero A (LCDM), and 3) open CDM model (OCDM). Using 
these samples, we show that the mass functions of rich clusters are scaling invariant. In §4, the 
predictions of scaling are compared with observations of clusters identified from IRAS sample. In 
§5, the index of the scaling is detected from samples of X-ray and gravitational lensing clusters 
with mass estimates. Finally, we describe our main findings in §6. 
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2. Theoretical background 

It is generally believed that cosmic gravitational clustering after the decoupling can roughly 
be described by three regimes: 1. linear regime; 2. quasilinear regime which is dominated by 
scale-invariant radial infall; and 3. fully developed nonlinear regime dominated by virialized 
nonradial motion. Since the amplitudes of the cosmic temperature fluctTiations revealed by COBE 
arc as small as AT/T ~ 10~^ (Bennett ct al. 1996), the gravitational clustering should remain in 
the linear regime, or at most in the quasilinear regime, on comoving scales larger than about 10 
h^^ Mpc and at redshifts higher than 2. On the other hand, observations indicate that clusters 
of galaxies with size of about 1.5 h~^ Mpc are probably the largest fully developed or virialized 
objects. Thus, the quasilinear evolution should be substantial on the scales from about 2 to 10 h~^ 
Mpc. Obviously, these facts are useful for discriminating among models of structure formation. 
Because different regime generally has different behavior of scaling, an effective way of doing 
model discrimination would be to compare the model predicted scaling with observations. 

In linear regime, the scaling is straightforward. If the primeval density perturbations are 
Gaussian and have power law spectrum P{k) = \5{k)\'^ = Ak'^\ the variance of the initial mass 
fluctuation within a spherical region on length scale R is given by {M)ji oc i?(3-")/2. Therefore, 
there is a scaling of i? ^ R', and {M)n, {R' / R)^^-''^^ (M) r. Since {M)r is a statistical average 
with respect to N{M, R)dM, the number density of radius R halo with mass from M to M + dM, 
the scaling of {M)ji indicates the scaling of the number density as 

n[> (i?7i?)-(^-")/2M, R'] = n{> M, R) (1) 

where 

n{>M,R)= N{M,R)dM (2) 
Jm 

which is the accumulative number density of radius R halo with mass larger than M. 

Since R is the initial radius of the halos, not the radius rd of the developed halos, the scaling 
Eq.(l) cannot be directly tested by observed halos. Only in the stage of linear regime, the scale of 
halos rd oc R, we have then scaling relationship as 

n[> ir'd/rd)-''M, r'^] = n(> M, r^) (3) 

where scaling index = (3 — n)/2, subscript L is for linear regime. Namely, the scaling index a 
is a function of n, or it is model-dependent. 

Beyond linear regime, the existence of scaling behavior of gravitational clustering can be 
shown by comparing the gravitational clustering of cosmic matter with the evolution of the 
profile of a growing interface, or surface roughing. The later is described by the so-called KPZ 
equation (Kardar, Parisi &; Zhang 1986), which mainly consists of terms of nonlinear evolution and 
stochastic force, or stochastic initial perturbations. A major breakthrough in this approach has 
been to show the existence of universal scaling properties of the surface growth by the dynamical 
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renormalization group (Barabasi & Stanley 1995)). It has been found that the equations describing 

the evolution of cosmic density perturbations are essentially KPZ-like (Berera & Fang 1994). The 
equations of cosmic gravitational clustering contains the similar nonlinear term and a stochastic 
initial perturbations. The gravitational potential of cosmic matter plays the similar role as the 
height h{x) of the surface. Thus, the structure formation in the universe substantially is also 
a phenomena of structural "surface" growth: an initially flat or smooth 3-dimension surface 
described by the Robertson- Walker metric evolved into a roughen one. Consequently, the cosmic 
gravitational clustering should also show the scaling feature as the KPZ-like surface growths. 

The scaling of a KPZ-like surface is given by ((/t^ - h'^)) oc L^f{t/L^), where L is the length 
of the coarse-grained average; t is time. The function f{x) is ~ x^^^ when x is small, and f{x) ^ 
const, when x is large (Vicsek 1992). The scaling indexes x z depend on the indexes of the 
spectrum of the stochastic force or the initial perturbations on the surface. Because Vd is, in fast, 
the scale of a coarse-grained average of mass field of the universe, we have 

(M)«r3/(t/rf,) (4) 

Since the scaling indexes a and /? depend on the indexes of the power spectrum of the initial 
perturbations, they are model-dependent. 

From eq.(4) one can generally conclude the existence of scaling during quasilinear and 
non-linear regimes. The scaling is characterized by 1.) for small scale r^/ and/or longer (later) time 
t, the coarse-grained mass distributions have scaling (M) oc r^; 2.) for larger r^i and/or earlier 
time t, the mass distribution will deviate from the a scaling. Because larger and earlier time 
t corresponds to quasilinear regime, the deviation from a-scaling should be due to the difference 
between the scaling behaviors of the quasilinear and non-linear regimes. The indexes a and /3 are 
initial-perturbation dependent. Thus, one can expect that the a-scaling and the deviation from 
this scaling are useful for discriminating among models of structure formation. 

These general conclusions are illustrated by some semi-analytical approaches. For instance, 
with the assumption of "stable clustering", the index 7 of the two-point correlation function 
.^(r) oc r~'^ at the fully developed non-linear regime is found to be 7 = 3(n -|- 3)/(n -|- 5) (Peebles 
1965). This turns to a = 3/(n -|- 5) (Pandmanabhan h Engineer 1998). This means that the 
a-scaling of highly virialized halos is indeed dependent on the initial density perturbations. 

As for the deviation from the a-scaling, one can refer to the power transfer via mode coupling 
in quasilinear regime. It has been found that for CDM-like spectrum, the power transfer of density 
perturbations is from large scales to small ones (e.g. Suto & Sasaki 1991). This is, the larger 
scale perturbations relatively have higher power in the quasilinear regime than fully developed 
non-linear regime. Therefore, the index a in quasilinear regime will generally be larger than that 
of non-linear regime. Thus, one can predict that the index a should show a "going up" on larger 
scales. This conclusion can also be illustrated by spherical in-fall model. Using scale-invariant 
spherical in-fall approximation, it has been shown that an initial scaling of may lead to a 
scaling of a' = 3/(4 — a^). For a = (3 — rz)/2, we have ol = 6/(n -|- 5), which is larger than the a 



-5- 



given by "stable clustering" . Therefore, the scaling index in the quasilinear regime is larger than 
that of non-linear regime (Padmanabhan 1996). 

Despite these semi-analytical results are consistent with the scenario of KPZ-like dynamical 
scaling in general, it is difficult to semi-analytically calculate the accurate relation between the 
scaling index and various initial perturbations. One cannot test models by the scaling index 
calculated from assumptions like the "stable clustering" or spherical in-fall. The studies of 
surface growth has shown that in the case of 3-D one cannot find the index of dynamical scaling 
analytically (Vicsek 1992). To find the number of the scaling, and to test models by this scaling, 
numerical study is necessary. This motivated us to investigate the dynamical scaling of the 
gravitational clustering numerically. 

3. Mass functions of clusters and its scaling 

3.1. Multi-scale identifications of clusters 

In order to study the mass functions of clusters at different scales, we have performed N-body 
simulations with the P^M code (Jing & Fang 1994; Jing et al. 1995; Wu et al. 1997b) for 
models of the SCDM, LCDM and OCDM. The cosmological parameters {i^M,^A,h,as) are taken 
to be (1.0,0.0,0.5,0.62), (0.3,0.7,0.75,1.0), (0.3,0.0,0.75,1.0) for the SCDM, LCDM and OCDM, 
respectively. The models with these parameters seem to provide a good approximation to many 
observational properties of the universe, especially the abundance of clusters (e.g. Jing & Fang 
1994; Bahcall, Fan & Cen 1997). 

The simulation parameters, including box size L, number of particles Np and the effective 
force resolution ij, are chosen to be {L,Np,ri) = (310/i~^ Mpc, 64:^,0.24h~^ Mpc). We have run 8 
realizations for each model. A particle has mass of 3.14 x 10^^Q,Mh~^MQ, which is small enough 
to resolve reliably the rich clusters with M > 3.0 x IQ-*^^^"^ Mq. 

To effectively identify clusters with different comoving radius, we haven't employed the 
traditional friends-of-friends algorithm, but instead developed an algorithm based on the discrete 
wavelet transformation (DWT) (see §Al of Appendix). The details of the DWT identification 
of clusters will be reported separately (Xu, Fang &: Deng, 1998). Briefly, wc first describe the 
distribution of the particles by a 3-D matrix, and then do fast 3-D Daubechies 4 DWT and the 
reversed transformations to find the wavelet function coefficients (WFCs) and scaling function 
coefficients (SFCs) on various scales. For each scale, the cells with SFCs larger than those of the 
random sample by a given statistical significance are picked up as cluster candidates. Around each 
of the candidates, we place a 6^ grid with the size of cluster diameter and search for the accurate 
center. The cluster center is taken as the position with largest mass surrounded. The mass is 
measured by counting the particles within a sphere of radius rd, the volume of which is the same 
as the cells. Whenever two clusters are closer than 2rci, the cluster with smaller M is deleted from 
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the list. We iterate the above steps for particles which have not been listed as cluster members 
until no further clusters are found. We will call clusters identified by radius rd as rcrdusters. 
Since the DWT technique treats tlic identification at different r^i in a uniform way, it is suitable 
to study the rcz-dependencc of clusters. 

Fig. 1 shows the derived cluster mass function n(> M, Tc/), which is the number density of 
clusters with mass larger than M within radius rd- For all the three models the mass functions of 
Tci = 1.5h~^Mpc clusters given by the DWT method are found to be in good agreement with those 
derived from friends-of-friends identification (Jing Sz Fang 1994). At higher z, these two methods, 
the friends-of-friends and DWT, also provide the same mass functions for Vd = 1.5^~^Mpc 
clusters. It turns out that the DWT technique of identifying clusters is indeed reliable. This 
algorithm can also be applied to real samples of galaxy distribution. Since the number densities 
of galaxies in real samples are much less than those of particles of N-body simulation samples, 
sampling error may lead to some false identification of poor clusters. However, for rich clusters, 
the effect of sampling error is small. 

Since in our identified sample each cluster is characterized by two parameters M and rd, it 
is inconvenient to define richness by mass alone. Instead, the relative richness can be defined by 
number density. This is, clusters with mass Mi and scale ri are considered to have the same 
richness as clusters with M2 and r2 if n(> Mi,ri) = n(> M2,r2). In this paper, the number 
density n is expressed as where d is the mean separation of the clusters. Therefore, clusters 

with n{> M,rd) < l(50/i~^ Mpc)~^ correspond to rich clusters of M > 5.5 x 10^^/i~^ Mq on scale 
rd = 1.5/1-1 Mpc for SCDM at z=0. 



3.2. The scaling of mass functions 

It can be easily seen from Fig. 1 that all the mass functions on various scales have a similar 
shape. This is, if the 1.5 h"^ Mpc mass functions are shifted horizontally along M-axis, they can 
approximately coincide with the mass functions of rd = 3, 6 and 12 Mpc clusters, respectively. 
Especially, a good match is found in the range of abundances lower than, or richness higher than, 
1 (50 h-^ Mpc)-3. 

This similarity can indeed be described by a scale invariance of the abundance mentioned in 

§2 

n[> {rd/roTM, ra] ^ n{> M, ro) (5) 

where a is the index of the scaling. Eq.(5) indicates that the number density of halos (clusters) 
having masses > M within radius ro is the same as that of halos having masses > {rd/ro)'^M 
within rd- We have tested eq.(5) by calculating a for n(> M,rd) in the abundance range of 
n{> M,rd) < 1(50/1-1 Mpc)-^, and found that a remains roughly to be a constant within 5 — 10% 
in the radius range rd ^ I — 6h~^ Mpc when ro = 1.5h~^ Mpc. 
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The scale-invariance is conveniently expressed by a mass-radius scaling which is the solution 
of the following equation: 

n[> M{r,i), r,i] = n[> M(1.5), 1.5] (6) 
where 1.5 denotes ra = 1.5h~^ Mpc. If eq.(5) is correct, we will have 

M(r,0 = (^)"M(1.5) (7) 

(In this paper, we denote M from scaling by M{rci), the mass within r by M(r).) 

In Fig. 2, we plot log[M(rci)/M(1.5)] against log from the solution of eq.(6) for the three 
models, in which the "richness" is taken to be n[> M(1.5),1.5] = l{90h^^ Mpc)^"^. The loci 
of log[M(rc;)/M(1.5)] vs. logr^i can be fairly well approximated by straight lines in the range 
of 1.5 < rci < Qh~^ Mpc. Especially, for the LCDM the relation of log[M(rci)/M(1.5)] - logra 
is quite straight in this range. For SCDM and OCDM they also do not deviate too much from 
straight lines. Therefore, the slopes of these lines, or a's, nearly remain constant in the radius 
range 1.5 < r^i < Qh"^ Mpc and "richness" range n[> M{rci),rci] < l{50hr^ Mpc)~^. The straight 
line fitting yields a «i 0.60 for LCDM and OCDM, while a 0.80 for SCDM. 

Fig. 2 shows that all curves log[M(rci)/M(1.5)]-log are going up around rd ~ 6h~^ Mpc. 
In other words, the index a becomes larger in the range of > 6h~^ Mpc. This is qualitatively 
in good agreement with the behavior of the scaling predicted in last section. As the gravitation 
clustering on scales larger than 6h~^ Mpc probably still partially remains of infall evolution, the 
scaling index is affected by the quasilinear evolution, and higher than that on scales less than 
6^~"^Mpc. Indeed, on such large scales, more and more (M, rc;)-identified halos don't have dense 
virialized cores in their centers, but arc with complicated and irregular shapes, one of which is 
illustrated in the Fig. 1 of Jing and Fang (1994). We also find that the "going up" behavior of the 
scaling index a is more significant at the quasilinear regime at earlier times. A similar result can 
also be seen in the comparison of universal profile with N-body simulation on larger scale (Diaferio 
& Geller, 1997). 



3.3. The "universal" density profile 

It has been proposed that the mass density profiles of regular clusters may universally be 
fitted by M(r) oc ln(l + r/a) — r/(r + a), where a is proportional to the virialized radius of clusters. 
Both r and a arc in unit of r2oo, the radius where the average interior density is 200 times higher 
than the critical density of the universe (Navarro, Frenk & White 1996). In the range of r > a, 
the universal mass profile can be approximated as a power law 

M(r) (X r"^'-" (8) 

Obviously, the index apro depends on the fitting of developed virialized core with a (Carlberg et 
al. 1997b). 
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Eq.(8) looks very similar to eq.(7). However, we should not identify the index apro with a, 
because in eq.(7) is the scale used to identify clusters, while r in eq.(8) is the radius of each 
cluster. Opro will be the same as a only if clusters identified by different rd have the same mass 
density profile. Fig. 3 shows, however, a systematic difference of mass density profile of clusters 
with different rd- Moreover, the larger the scale, the larger the systematic difference. This is 
because the clusters identified by large rd are still partially in the quasilinear regime for which the 
halos are pre-virialized. Halos, which lack a virialized core in their centers, cannot be described 
by the universal density profile. 

It should be pointed out that an original meaning of the "universal" is the independence of 
profile on the initial condition, universal number regardless the parameters of initial 

spectrum like the index n. Namely, all initial parameters are forgotten in the late time (non-linear 
regime) clustering. Theoretically, it is far from clear of the condition under which the late time 
profiles do not remember the initial condition. At least, the assumption of the stable clustering 
cannot co-exist with an initial-condition- independent profile (Pandmanabhan & Engineer 1998). 
Fig. 3 shows that the mean profiles M(r) of clusters are significantly different for different models. 
Namely, the late time profiles do remember the initial condition. 

4. Scaling invariance of IRAS clusters 

4.1. Multi-scale identified clusters from IRAS 

We can directly test the scale invariance via eqs.(6) and (7). To do this, we have applied 
the multi-scale DWT method to identify galaxy clusters in the redshift survey of the Infrared 
Astronomical Satellite (IRAS) galaxies with flux limit of fao > 1.2 Jy (Fisher ct al. 1995b). The 
IRAS surveys arc uniform and complete down to galactic latitudes ±5° from the galactic plane. 
The 1.2 Jy IRAS survey consists of 5313 galaxies, which cover 87.6% of sky, and with 12.34% of 
the sky belonging to the so-called "excluded zones" . 

We fill up the "excluded zones" and holes with a randomly distributed galaxies, which are 
generated with the same mean number density, and the same radial (redshift) selection function as 
other areas. This treatment may lead to a little underestimate of the number density of clusters. 
However, in terms of the ratio between abundances of clusters on different scales the effect of 
"excluded zones" is small. Moreover, random data, like Poisson distributions, may give rise to 
false statistics on order higher than second, but this problem will not be met in the statistics of 
abundance. 

As usual, to reduce the effect of radial selection function, we divide the sample into a series of 
shells with thickness ~ 2500 km s~^ in redshift space, and measure the number count functions 
n(> Ng,rd) in three shells of v: [2500-5000], [5000-7500] and [7500-10000] km s"^ There are only 
few galaxies beyond 10000 km s~^. They are not suitable for studying the scaling of the cluster 
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abundance. 

With these redshift shells Zi the density field 5p{s, Zi) can be reconstructed by a 2-D DWT 
decomposition with bases V'j,i(s) as 

(5p(s, Zi) = Y,Y. *j,iV'j,i(s) (9) 

3 1=0 

where the WFCs w^^i can be found from the galaxy distribution of the sample (A2 Appendix). 

Since the wavelet bases ■i/'j,! are orthogonal and complete, the density field 5p{s, Zi) is the same 
as descriptions by other orthogonal-complete decomposition, say the bases of spherical harmonics 
and Bessel functions (Fisher et al. 1995a). It has been shown that the power spectrum detected 
by the DWT is the same as that by a Fourier decomposition (Pando & Fang 1998). The mass 
density given by the DWT is also the same as that of Bessel-spherical harmonics. We choose the 
DWT bases only because it is easy to detect the j-dependence, or scale-dependence. 

We will use 5p{s, Zi) to identify clusters on various scales j by the same way developed in 
§3.1. Similar to the N-body simulation, the cluster center is defined as the point around which a 
maximum number of galaxies arc enclosed within a cylinder with length of 3000 km s~^ and radius 
equal to that of the cluster. Whenever two clusters are closer than 2rcu the cluster with smaller 
number of galaxies Ng is deleted from the list. We iterate the steps until no further clusters are 
found. The result gives a number count function njjiAs{> ^g-,^cl)-, which is the number density of 
IRAS clusters consisting of > Ng galaxies within radius rd- 



4.2. Ng-M conversion 

Mass is not directly measurable for IRAS clusters. To transfer the number count function 
n{> Ng,rci) to mass function n(> M,rci), we need a conversion between the number count and 
mass for the IRAS clusters. Because the total luminosity of a ri.5-cluster is proportional to its 
richness (mass), and the mass-luminosity ratio for these cluster is independent of the total mass 
of the clusters (Bahcall & Cen 1993), it is reasonable to assume that the mass M of a cluster is 
proportional to the total number Nfotal of the galaxies in the given cluster, i.e. 

M = ANtotal (10) 

Obviously, A is only a number-mass conversion coefficient, not the mean mass of galaxies, as 
the cluster mass M is dominated by dark matter. The total number of galaxies is related to the 
counted galaxy by the selection function ^(2;) (Fisher et al. 1995b): 

Ntotal = Ng/4>{z) (11) 



We have then 



M = A* Ng/4){z) 



(12) 
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We calibrate the coefficient A at a fixed mass, i.e., at a given abundance n 10 ^h^ Mpc ^ 
tfirougfi tfie equation: 

niRAs{> Ng, 1.5) = nBc{> M, 1.5) (13) 

wfiere subscript BC means the mass function of Bahcall and Cen (1992), and niRAs{> Ng, 1-5) is 
from the shell of [5000-7500] kms~^. We solve for Ng and M. Using this pair of Ng and M, we 
derive the coefficient A from eq.(12). 

With this A, one can transfer the number count function niRAs{> Ng,Tcl) into mass function 
nijiAs{> M,rci) for the entire range of Ng. The result of njjiAs{> M,rci) is plotted in Fig. 4, in 
which the horizontal errors of IRAS clusters are caused by the Poisson errors of number counting 
of galaxies in a given cluster, and the vertical errors are from the Poisson errors of counting the 
clusters. 

The BC mass function is also shown in Fig.4. It has been known that SCDM model doesn't fit 
EC's MF while LCDM fits well(Jing & Fang 1994). IRAS results confirm this conclusion. One can 
find that the mass function of the IRAS clusters with rd = 1.5/i~^ Mpc is basically the same as 
the mass function of BC sample, especially, for clusters with richness M > Mq. Namely, 

the masses of clusters identified from the 1.2 Jy IRAS samples are statistically the same as those 
of the clusters in the BC sample. This result is consistent with the following fact: the IRAS 
galaxies trace the local large scale structures. It has been found that many optically identified 
structures, including super clusters, voids, great attractor and Abell clusters, have been identified 
from density field constructed from IRAS data (Webster, Lahav Sz Fisher 1997). Considering 
that the clusters of BC sample consist of optical and X-ray clusters, and the IRAS galaxies are 
biased, containing more later type galaxies, Fig. 4 indicates that the early-type galaxies map 
about the same mass field as late-type, despite the early-type galaxies are clustered more strongly 
than late-type galaxies. This is because that in terms of second order of statistics the segregation 
between the early and late-type galaxies is almost linear, and scale-independent at relatively large 
scales. Using Stromlo-APM redshift survey it has been shown that for second order statistics the 
scale-dependence of the segregation between early and late-type galaxies is not large than 1 cr in 
the range from 1 to 20 h'^ Mpc (Lovcday et al. 1995, Fang, Deng & Xia 1998). 

One can further test the reasonableness of the calibration (12) by comparing different redshift 
shells. With eq. (12) and the selection function 4>{z) (Fisher et al. 1995b), we found that the best 
values of A for the three shells are, respectively, 10^2•^=^°■^ io^^-^=^°-2, IQ^^ '^^o.s j^-i y^^_ r^hey 
are indeed the same within error limit. Thus, the mass functions niRAs{> M,rci) from the shells 
of [2500-5000] and [7500-10000] kms"^ are the same as [5000-7500] km s-\ and also the same 
as the BC sample. Therefore, the number of A provides a consistent Ng — M conversion for the 
entire sample of the 1.2 Jy IRAS galaxies. Using the conversion of Abell richness to cluster mass: 
M/Nc = 0.6 X 10i3/i-i Mq (Bahcah & Cen 1993), we have A = M/Ntotal = 0.8 x lO^^/j-i M© 
gives Nc/Ntotal ~ 7. It means every count of the 1.2 Jy IRAS galaxies (after selection-function 
correction) corresponds to about 7 times Abell count of optical galaxies in counting the mass of a 
cluster. 
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4.3. Scaling of IRAS clusters and models 

Because M/L reaches a constant asymptotic value beyond ^ 1 Mpc (Bahcall, Lubin k. 
Dorman 1995) , and there is no evidence of significant scale dependence of bias factor of IRAS 
galaxies from 1 to 10 h~^ Mpc, the Ng-M conversion eq.(12) should be applicable on scales 
> 1.5/i-^Mpc. Thus we can find mass functions n(> M,rci) from number-count functions 
n(> Ng,rci) of IRAS cluster with rd > 1.5/i"^Mpc. 

Using these IRAS n(> M,rci), we test the scaling by eq.(6). The solutions of eq.(6) for both 
IRAS data and simulation sample are shown in Fig. 5. The theoretical curves in Fig. 5 are similar 
to Fig. 2, and the richness parameter is taken to be l(50/i~^ Mpc)~^. 

Since only mass ratios M(rc/)/M(1.5) of the solutions eq.(6) are plotted in Fig. 5, the result 
doesn't depend on the value of A. The effect of linear bias of galaxies will also be canceled in 
this ratio. The errors of IRAS data in Fig. 5 are calculated from both horizontal and vertical 
errors of the mass function (Fig. 4). Since the mass function is very steep for rich clusters, i.e. 
\d\nn/d\nM\ S> 1, an uncertainty of Inn will transfer to a relatively small uncertainty of InM. 
Because the mean number density of 1.2 Jy IRAS galaxies is low, the major source of errors in 
Fig. 5 is from Poisson errors of the number of galaxies in clusters. 

Despite the errors are large. Fig. 5 already shows that the IRAS data is basically consistent 
with the prediction: there is a scaling invariance in the range of 1.5 < r^i < 4.5/i~^Mpc with index 
a ~ 0.5, and the scaling index of scaling is "going up" on larger scales. The results of simulation 
samples show that the index of scaling is model-dependent. The three panels of Fig. 5 generally 
are in agreement with the two low mass models (LCDM and OCDM) within 1 a, but show a 
systematic disfavor of the SCDM at > la. The error bars in Fig. 5 are slightly overestimated by 
assuming that the Poisson errors of mass estimates between two scales are independent. They 
might be about l/\/2 times smaller if the two Poisson errors are completely correlated (cloud 
in cloud scenario). So, in each of the three shells the SCDM is away from the data at the level 
of > 1 ~ 1.4(7, or disfavored at a confidence level of > 68% — 82%. If all the shells are binned 
together, the confidence should certainly be higher than that of individual shell, because the 
Poisson errors will be smaller. However, it is difficult to estimate the infiuence of the selection 
function. Moreover, if spiral galaxies were underrepresented within about 1.5 Mpc, the 
true values of the scaling index should be lower than that shown in Fig. 5, i.e. it strengthens 
the conclusion of disfavoring the SCDM. Therefore, the number 68% — 82% can be applied as a 
underestimated confidence level. 



5. Other detections of the scaling index 
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5.1. Sample of rich clusters with mass estimation 

Similar to eq.(2), the number density of clusters with mass in the range of M to M + dM 

is 

N{M,rci)dM = -^n{> M,rd)dM. (14) 
One can then define a mean mass of rd clusters with richness n(> M, rd) < uq by 

M{rd)= I MN{M,rd)dM. (15) 

Jn{>M,rci)<no 

Using eqs.(5), (6) and (7), we have 

log M{rci) = a log ra + const. (16) 

This means, the index a determined by fitting cq.(16) with a sample of mass measurements of 
clusters with richness n(> M,rci) < no should be the same as that given by abundance solution of 
eqs.(6) and (7). This prediction can be tested if we have fair samples of masses of clusters with 
richness n(> M, Vd) < uq and over a given radius range of Vd- 

Despite we still lack cluster mass samples covering a large radius range of rd and with the 
desired completeness, the current data are already able to preliminarily test the prediction eq.(16). 
For instance, it is generally believed that the weak gravitational lensing clusters being studied 
are among the richest clusters. Weak lensing mass estimate gives only a lower bound to the total 
cluster mass because of the unknown mean density in the so-called control annulus (Fahlman et al 
1994). Nevertheless, the radius dependence of the cluster mass given by weak gravitational lensing, 
Myji{rd), is insensitive to the control field which contributes only a constant component to the 
mass distribution. Moreover, under the isothermal assumption, we have M{rd)/rd on af, where 
(Ji is the sight-of-line velocity dispersion of the cluster galaxies. Since the relation M{rd)/rd oc erf 
is independent of richness, a velocity dispersion normalized mass, Mwi{rd)/cr1, appears to be less 
dependent on richness. So, the cluster mass sample consisting of Myji{rd) and velocity dispersion 
measurements would be good for fitting with eq.(16). There are 9 clusters which have both 
velocity dispersion measurements and weak gravitational lensing mass M^i estimates in the radius 
range of 0.15 to 2 hT^yipc (Wu & Fang 1997a). Although My,i{rd) from the weak lensing actually 
corresponds to a projected mass within rd-, the values of a given by either the projection or the 
3-D masses will be roughly the same at large radius Vd > lh~^ Mpc. A best fitting of this weak 
lensing data yields a = 0.71 ± 0.20. 

Another data set of rich cluster masses can be selected from X-ray measurements. The largest 
sample of X-ray clusters with mass estimates published to date is given by White, Jones & Forman 
(1998), which contains 226 X-ray cluster masses for 207 clusters derived from a deprojection of 
Einstein Observatory X-ray imaging data. Meanwhile, by extensively searching literature there 
are additional 152 X-ray determined cluster/group masses available. These data were obtained by 
either the similar approach to WJF (e.g. White &; Fabian 1995; Ettori, Fabian &; White 1998; 
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etc.) or the analysis of the ROSAT PSPC X-ray observations (e.g. Pildis. et al 1995; David, Jones 
& Forman 1995; Cirimele et al. 1997; etc.). All the available 144 measurements of cluster masses 
from X-ray observations with > 0.5/i~^Mpc are plotted in Fig. 6. 

The X-ray and optical measurements of some gravitational lensing clusters have shown that 
the gravitational lensing clusters on average correspond to X-ray temperature T > 7.5 keV and 
velocity dispersion of ai > 1200 km s^^. Therefore, it is reasonable to select richest clusters by 
the conditions of T > 7 keV and ai > 1000 km s~^. There are 11 measurements satisfying these 
conditions. These data are plotted as star-added circles in Fig. 6. It yields a = 0.66 it 0.40 in the 
range of 1.5 < ra < 4:h~^Mpc. 

It is interesting to see from Fig. 6 that the X-ray mass distribution has a clear upper envelope, 
and all the clusters selected by T > 7 keV and ai > 1000 km distribute along the envelope. 
It implies that the clusters at the envelope are among the richest clusters at the given r^/. Recall 
that observations may easily lose less massive clusters, but tend to pick up the bright and massive 
ones. This selection effect is severe for producing a reliable sample of less rich clusters, but benefit 
to the completeness of sample of rich clusters. Therefore, one can employ the envelope clusters 
to fit with eq.(16). To this end, we have binned the data set of cluster masses into 8 logarithmic 
intervals according to radius from Vd = 0.5 — 4.0/i~^Mpc. Within each bin, the mean value of the 
second and the third largest clusters is used as the envelope value. Two of the 8 bins having less 
than 2 clusters are not considered as reliable measurements of the envelope and are thus not used. 
All envelope values are shown in Fig. 6, in which the vertical error bars are the mass difference 
between the second and third masses of the most massive clusters. The least-square fitting gives 
a = 0.52 lb 0.25 in the radius range of 0.5 to 3 h~^Mpc. To reduce the possible effect of the bin 
size selection, we re-calculated a using different bin sizes. It turns out that differences among 
these results are not larger than la. 

All the detections of eq.(16) with independent ensembles of rich cluster mass estimates 
have yielded the same number of a within error limits. For clusters with "richness" 
n(> M,rci) < l{50hr^ Mpc)~^ in the range of 1 < r^/ < 4h~^ Mpc the mean value of a is 

0.63 lb 0.10. These values of a basically agree with the a detected by the abundance of IRAS 
clusters (§4.3). Therefore, current data are in good consistence with the scaling invariance of rich 
cluster abundance. 



5.2. Deviation from a constant a 

As has been discussed in §2, the scaling index a remains as a constant on smaller scales and 
later time, but will show "going up" on larger scales and early time on which the quasilinear 
evolution still plays important role. Therefore, one can expect that the "going up" behavior will 
be more significant at higher redshifts. An effective measure of the "going up" is the abundance 
ratio n[> {rci/ro)'^M,rci\/n{> M, ro) at higher redshifts. If the scaling invariance of abundance 
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perfectly holds with a constant index a, we have 



n{> M,ro) 



= 1 



(17) 



Therefore, the radius r beyond which this ratio no longer remains equal to unity is a measure of 
the importance of the quasilinear evolution. 

Fig. 7 plots the ratio n[> {rci/ro)°'M,rci]/n{> M,ro) of simulated samples with parameters 
a = 0.63, ro = 1.5/i~^ Mpc and M = 5.5 x lO^^/i^^M©. Fig. 7 shows that all the abundance ratios 
become larger than 1 on large scales. Namely the real values of a for the simulated samples are 
larger than the assumed 0.63 on larger scales. This is the "going up", indicating the deviation 
from constant a. 

Fig. 7 shows that this "going up" behavior is significantly diff'erent for different models. The 
SCDM curves are much more quickly and strongly "going up" than LCDM and OCDM. This is 
because clusters in the SCDM formed later than in the LCDM and OCDM. In the latter case, the 
ratio n[> {rci/ro)"M, rci]/n{> M, ro) remains equal to about unity in the range of rg = 1.5 — 6h~^ 
Mpc and z < 0.8, while for the former the corresponding radius range is much smaller. We have 
re-calculated Fig. 7 by changing the value of a. The results show that in the range a ~ 0.43 — 0.77, 
the OCDM and LCDM always have a larger radius range of the invariance than SCDM. Therefore, 
the radius range in which the abundance scaling invariance holds is effective for discriminating 
models. 

The fitting done in §5.1 has shown that a constant a of 0.52 - 0.70 in the radius rg from 1 to 
4 Mpc is consistent with all data. Moreover, considering most gravitational lensing and X-ray 
clusters used in §5.1 are of moderate redshift, the numbers of a 0.52 — 0.70 and 1 < rg < 4h~^ 
Mpc are actually real for moderate redshift. Thus, the SCDM will be in difficulty with the 
measurements of scaling invariance. It should be pointed out that the radius range 1 < rd < 4hr^ 
Mpc in §5.1 is for the physical scale of relevant clusters. A physical radius range 1-4 Mpc 
corresponds to a comoving range of 1.5 - 6 Mpc at redshift ~ 0.5. Therefore, this test is more 
robust at higher redshift. 

It may be difficult to directly test the scaling at high redshifts, because mass determination 
is often limited by the luminosity(or surface brightness) detection threshold of the survey. In 
this case we may turn to the correlation function. The scaling gives testable prediction on the 
correlation lengths of clusters on different scales. This will be investigated in Xu, Fang &; Deng 



(1998). 



6. 



Conclusions 



We have studied the scaling invariance of abundance of rich clusters. Both the N-body 
simulation and the available observational data have confirmed the existence of the scaling 
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invariance. This scaling is characterized by index a. which can be determined by the abundances 
of clusters on different scale r^;, or by fitting mass-radius relation of clusters. The scaling gains a 
further support from the following result: the a given by X-ray and gravitational lensing mass 
estimates is the same as that from the IRAS cluster abundance. Despite the significance level of 
current results is not high, it is worth of revealing the general behavior of the scaling of cosmic 
clustering in different evolutionary stages. It can already be employed for discrimination among 
models. 

The LCDM and OCDM abundances can always have a scaling in a larger scale range 
(~ 1 - 6/i~^ Mpc) with a ~ 0.5 — 0.7 for z < 0.8. While the SCDM can only provide a smaller 
scaling range (~ 1 - 3/1"^ Mpc) with a <~ 0.8 and z < 0.8. If a is allowed to be > 0.80, the 
SCDM can provide a scaling in the range of 1 - 4 h^^ Mpc for z < 0.8. However, the SCDM is 
difficult to fit a ~ 0.50 — 0.70 and ~ 1 — 4/i~^ Mpc simultaneously. This result may cause some 
troubles for the SCDM. 

In general, the mass density profiles of clusters in the low density models are steeper than that 
of the corresponding clusters in higher Qm models (Jing et al. 1995). In other words, the index a 
for the CIm = 1 model is always larger than that in a low density model {CIm < 1)- Therefore, the 
current result of the mass-radius scaling seems to favor models with a lower mass density. 

WX thanks the World Laboratory for a scholarship. XPW was supported by the National 
Science Foundation of China and the National Science Council of Taiwan. 



Let us briefly introduce the DWT analysis of large scale structures, for the details referring to 
(Fang &; Pando 1997; Pando & Fang 1996, 1998, and references therein). We consider here a 1-D 
mass density distribution p{x) or contrast d{x) = [p{x) — p]/p, which are mathematically random 
fields over a spatial range < x < L. It is not difficult to extend all results developed in this 
section into 2-D and 3-D because the DWT bases for higher dimension can be constructed by a 
direct product of 1-D bases. 

Like the Fourier expansion of the field 5{x), the DWT expansion of the field S{x) is given by 



A. The DWT decomposition and reconstruction 



A.l. 



The discrete wavelet transform 



oo 2^-1 




(Al) 



j=0 1=0 



where tpj^i{x), j = 0,1,..., I = 0...2^ — 1 are the bases of the DWT. Because these bases are 
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orthogonal and complete, the wavelet function coefficient (WFC), wj^i, is computed by 

'^3,1 = j Six)'4}j,i{^)dx. (A2) 

The wavelet transform bases tpj^i^x) are generated from the basic wavelet 'il){x/L) by a 
dilation, 2-', and a translation /, i.e. 

'^j,l{^)=W\ ^{2^x/L-l). (A3) 

The basic wavelet ■0 is designed to be continuous, admissible and localized. Unlike the Fourier 

bases cxp{i2TTnx / L) , which are non-local in physical space, the wavelet bases i]jj^i{x) are localized 
in both physical space and Fourier (scale) space. In physical space, ipj^ilx) is centered at position 
lL/2^, and in Fourier space, it is centered at wavenumber 2it x 2^ /L. Therefore, the DWT 
decomposes the density fluctuating field S(x) into domains j, I in phase space, and for each basis 
the corresponding area in the phase space is as small as that allowed by the uncertainty principle. 
WFC Wj^i and its intensity \wj^i\'^ describe, respectively, the fluctuation of density and its power 
on scale at position lL/2^ (Pando & Fang 1998). 



A. 2. Reconstruction of density field 

Using the completeness of the DWT basis, one can reconstruct the original density field from 
the coefficient wj i. To achieve this, DWT analysis employs another set of functions consisting 
of the so-called scaling functions, (pj^i, which are generated from the basic scaling (p{x/L) by a 
dilation, 2-', and a translation Z, i.e. 

^jA^)=[j\ <P{2'x/L-l). (A4) 

The basic scaling is essentially a window function with width x/L = 1. Thus, the scaling 

functions (j)j^i{x) are also windows, but with width (l/2-')L, and centered at lL/2^. The scaling 
function (pj^i(x) are orthogonal with respect to the index I, but not for j. This is a common 
property of windows, which can be orthogonal in physical space, but not in Fourier space. 

For Daubechics wavelets, the basic wavelet and the basic scaling are related by recursive 
equations as (Daubechies 1992) 

<P{x/L) =Y.iai<l){2x/L-l) 

^{x/L) =Y.ihi^i2x/L + l) ^ ^ 



where coefficients ai and hi are different for different wavelet. In this paper, we use the Daubechies 
4 wavelet (D4), for which oq = (1 + V3)/4, ai = (3 + V3)/4, 02 = (3 - V3)/4, ag = (1 - V3)/4. 
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Prom Eq.(A5), one can show that the scahng functions (f>j,i{x) are always orthogonal to the 
wavelet bases ijjj'^i'{x) if j < j', i.e. 



j(pj^i{x)iljf^i>(x)dx = 0, for j<f. (A6) 

Therefore, (f)j,i{x) can be expressed by 'tpji^ii{x) as 

oo 2^' -1 j-1 2^' -1 

hA^) = II Cjl;3'l''4^j',l'{x) = Z Z Cji.fi>ijj',iix). (A7) 

j'=0 l'=0 j'=0 l'=0 

The coefficients Cji-j'i' = J (j)j^i{x)il>ji^i'{x)dx can be determined by ai and bi. 
Using (f)j^i{x), we construct a density field on scale j as 

23-1 

p>{x) = Wj,i(j)j,i{x), (A8) 
/=o 

where Wj^i is called the scaling function coefficient (SFC) given by 

'^hl = / p{x)(pji{x)dx. (A9) 
Jo 

Since the scaling function (t)j^i{x) is window- like, the coefficient Wj^i is actually a "count-in-cell" in 
a window on scale j at position I. 

Using Eqs.(Al), (A7), (A8) and (A9), one can find 

j'=0 l'=0 

Namely, p^{x) contains all terms of density fluctuations wy ^it{x) of j' < j, but not terms of 
j' > j- Prom Eqs.(Al) and (AlO), we have 

00 2J'-1 

p{x) = p>{x) + pJ2 Y Wf^i't/jf^i'ix). (All) 

j'=j l'=0 

One can also define the smoothed density contrasts on scale j to be 

S'ix) ^ ^^^4-^ = E (A12) 

P j'=0 l'=0 

Eqs.(All) and (A12) show that and are smoothed density fields on scale j. One can 

construct the density field p^^x) or 6^x) on finer and finer scales by WFCs wj^i till to the precision 
of the original field. Since the sets of bases tpj^i and (f)j^i are complete, the original field can be 
reconstructed without lost information. 
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Fig. 1. — Mass functions of clusters identified with radii = 1.5, 3.0, 6.0 and 12 h~^Mpc for 
SCDM, LCDM and OCDM models at 2; = 0. n(> M,rci) is the number density of clusters with 
mass larger than M within radius rd- M is in unit of h"^ Mq. The observed data are for clusters 
with radius = 1.5/i~^Mpc (BahcaU & Gen 1992). 

Fig. 2. — Mass-radius scaling of clusters given by the solution of abundance equation (6). The 
"richness" is taken to be n = 1(90^~^ Mpc)~^. The thin, dashed, and thick lines are for SGDM, 
LGDM and OGDM, respectively. 

Fig. 3. — The mean mass-radius relation M(rci) of clusters with different identification scales rd- 
The "richness" for all ra clusters is n = l{90h~^ Mpc)~^. The dashed and solid lines are for SCDM 
and OCDM, respectively. 



This preprint was prepared with the AAS IM^^X macros v4.0. 
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Fig. 4. — Mass function of rd = 1.5h ^Mpc IRAS clusters (circles with error bars) in the shell 
[5000,7500]km s'^. The stars are for the mass function of Bahcall & Cen (1992). The dashed line 
shows the mass function by FOFs method of Jing & Fang (1994). The thick line shows the mass 
function by DWT method averaged over 5-realizations of our 310/i~^Mpc box simulations. The thin 
line shows the mass function by DWT method from 1 realization of 155/t~^Mpc box simulation. 

Fig. 5. — Mass-radius scaling of clusters given by the solution of abundance equation (2). The 
thick, thin, and dashed lines are for LCDM, OCDM and SCDM, respectively. The data with error 
bars come from the clusters identified from IRAS galaxies. The richness parameter is taken to be 
l(50/i-i Mpc)-3. 

Fig. 6. — All the available 144 measurements (circles) of cluster masses from X-ray observations with 
rd > 0.5/i~^Mpc. The 11 richest clusters(r > 7 keV and a > lOOOkms"^) with ra > 1.5/i~^Mpc 
are marked additionally with stars. Squares with error bars are the data of upper envelope clusters. 
The horizontal error bars showing the widths of radius bins, and the vertical error bars being the 
difference between the masses of the second and the third most massive clusters within each bin. 
The solid line is an equal-weight least square fit to the envelope data. 

Fig. 7. — n[> Mth{rci),rci]/ni,5 vs. rd in the range of = 1.5 — 12/i""-'-Mpc at redshifts z ~ 0.2, 0.5 
and 0.8, for SCDM, LCDM and OCDM. Here Mth{rd) = {rd/n.^TMi,^ and ni.5 = n(> M1.5, 1.5). 
The parameters arc taken to be Mi. 5 = 5.5 x 10"*^^^"^ Mq, ri.5 = 1.5/i~^Mpc, and a = 0.63. The 
dot-dashed horizontal lines are for perfect scale invariance with a = 0.63. All curves for simulated 
samples show the "going up" of a on large scales. 
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